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ABSTRACT 

After a discussion of the properties of degenerate fermion balls, we analyze 
the orbits of the stars SO-1 and SO-2, which have the smallest projected distances 
to Sgr A*, in the supermassive black hole as well as in the fermion ball scenarios 
of the Galactic center. It is shown that both scenarios are consistent with the 
data, as measured during the last six years by Genzel et al. and Ghez et al. The 
free parameters of the projected orbit of a star are the unknown components 
of its velocity and distance z to Sgr A* in 1995.4, with the z-axis being in 
the line of sight. We show, in the case of SO-1 and SO-2, that the z — Vz phase- 
space, which fits the data, is much larger for the fermion ball than for the black 
hole scenario. Future measurements of the positions or radial velocities of SO-1 
and SO-2 could reduce this allowed phase-space and eventually rule out one of 
the currently acceptable scenarios. This may shed some light into the nature of 
the supermassive compact dark object, or dark matter in general at the center 
of our Galaxy. 

Subject headings: black hole physics-celestial mechanics, stellar dynamics dark 
matter - elementary particles - Galaxy: center 
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1. Introduction 

There is strong evidence for the existence of a supermassive compact dark object near 
the enigmatic radio source Sagittarius A* (Sgr A*) which is located at or close to the 
dynamical center of the Galaxy (Rogers et al. 1994; Genzel et al. 1997; Lo et al. 1998; 
Ghez et al. 1998). Stars observed in the 2.2 //m infrared K-band at projected distances ^ 
5 mpc from Sgr A*, and moving with projected velocities ^ 1400 km s~^, indicate that a 
mass of (2.6 ± 0.2) x IO^Mq must be concentrated within a radius ~ 15 mpc from Sgr A* 
(Haller et al. 1996; Eckart and Genzel 1996, 1997; Genzel & Townes 1987; Genzel et al. 
1994, 1996, 1999, 2000; Ghez et al. 1998, 2000). VLBA radio interferometry measurements 
at 7 mm wavelength constrain the size of the radio wave emitting region of Sgr A* to ^ 1 
AU in E-W direction and ~ 3.6 AU in N-S direction (Rogers et al. 1994, Bower and Backer 
1998, Krichbaum et al. 1994, Lo et al. 1998), and the proper motion of Sgr A* relative to 
the quasar background to ^ 20 km s~^ (Baker 1996; Reid et al. 1999; Baker and Sramek 
1999). As the fast moving stars of the central cluster interact gravitationally with Sgr A*, 
the proper motion of the radio source cannot remain as small as it is now for ~ 200 kyr 
unless Sgr A* is attached to some mass ^ lO^M©. In spite of these well-known stringent 
facts, the enigmatic radio source Sgr A*, as well as the supermassive compact dark object 
that is perhaps associated with it, are still two of the most challenging mysteries of modern 
astrophysics. 

It is currently believed that the enigmatic radio source Sgr A* coincides in position 
with a supermassive black hole (BH) of (2.6 ± 0.2) x lO^M© at the dynamical center of 
the Galaxy. Although standard thin accretion disk theory fails to explain the peculiar low 
luminosity ^ 10^'' erg s^^ of the Galactic center (Goldwurm et al. 1994), many models have 
been developed that describe the spectrum of Sgr A* fairly well, based on the assumption 
that it is a BH. The models proposed for the radio emission, range from quasi-spherical 
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inflows (Melia 1994; Narayan and Mahadevan 1995; Narayan et al. 1998; Mahadevan 1998) 
to a jet-like outflow (Falcke, Mannheim and Biermann 1993; Falcke and Biermann 1996; 
Falcke and Biermann 1999). Yet, as some of these models appear to contradict each other, 
not all of them can represent the whole truth. We also note that the Galactic center is a 
weak source of diffuse emission in the 2-10 keV energy range and in the lines of several ions 
(Sunyacv ct al. 1993; Koyama ct al. 1996; Sidoli and Mereghetti 1999). Thus, apart from 
earthbound VLBA radio interferometers, space missions such as the European Multi-Mirror 
satellite (XMM) and Chandra X-ray sateUite, may eventually provide conclusive evidence 
for the nature of Sgr A* and the supermassive compact dark object at the Galactic center. 
In fact, the Chandra X-ray satellite has recently detected a point source at the location of 
Sgr A* (Baganoff et al. 1999) with a luminosity two times smaller than the upper limit set 
by the ROSAT sateUite some years ago (Predehl and Triimper 1994). For more detailed 
recent reviews on the Galactic center wc refer to Morris and Serabyn 1996, Genzel and 
Eckart 1999, Kormendy and Ho 2000, and Yusef-Zadeh et al. 2000. 

Supermassive compact dark objects have also been inferred at the centers of many 
other galaxies, such as M87 (Ford et al. 1994; Harms et al. 1994; Macchetto et al. 1997) 
and NGC 4258 (GreenhiU et al. 1995; Myoshi et al. 1995). For recent reviews we refer to 
Richstone et al. 1998, Ho and Kormendy 2000, and Kormendy 2000. In fact, perhaps with 
the exception of dwarf galaxies, all galaxies may harbor such supermassive compact dark 
objects at their centers. However, only a small fraction of these show strong radio emission 
similar to that of the enigmatic radio source Sgr A* at the center of our Galaxy. For instance 
M31 does not have such a strong compact radio source, although the supermassive compact 
dark object at the center of M31 has a much larger mass (~ 3 x lO^M©) than that of our 
Galaxy (Dressier and Richstone 1988; Kormendy 1988). It seems, therefore, prudent not to 
take for granted that the enigmatic radio source Sgr A* and the supermassive compact dark 
object at the center of our Galaxy are necessarily one and the same object. 
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An unambiguous proof for the existence of a BH requires the observation of stars 
moving at relativistic velocities near the event horizon. However, in the case of our Galaxy, 
the stars SO-1 and SO-2, that are presumably closest to the suspected BH, reach projected 
velocities ^ 1400 km s~^. Assuming a radial velocity of — 0, this corresponds to the 
escape velocity at a distance <; 5 x 10^ Schwarzschild radii from the BH. Thus any dark 
object, having a mass ~ 2.6 x lO^M© and a radius ^ 5 x 10"^ Schwarzschild radii, would 
fit the current data on the proper motion of the stars of the central cluster as well as the 
BH scenario. One of the reasons why the BH scenario of Sgr A* is so popular, is that 
the only baryonic alternative to a BH that we can imagine, is a cluster of dark stars (e.g. 
brown dwarfs, old white dwarfs, neutron stars, etc.), having a total mass of ~ 2.6 x IO^Mq 
concentrated within a radius of ~ 15 mpc. However, such a star cluster would disintegrate 
through gravitational ejection of stars on a time scale ^ 100 Myr, which is much too 
short to explain why this object still seems to be around today ~ 10 Gyr after its likely 
formation together with the Galaxy (Sanders 1992; Haller et al. 1996; Maoz 1995, 1998). 
Nevertheless, in order to test the vahdity of the BH hypothesis meaningfully, we definitely 
need an alternative and consistent finite size model of the supermassive compact dark 
objects at the galactic centers. 

2. The case for degenerate fermion balls 

It is well known that our Galactic halo is dominated by dark matter, the bulk part of 
which must be nonbaryonic (Alcock 2000). Numerical simulations show that dark matter 
in the form of a gas of weakly interacting massive particles, will eventually produce a 
high-density spike at the center of the Galaxy (Navarro et al. 1997; Gondolo & Silk 1999). 
It is therefore conceivable that the supermassive compact dark object at the center of our 
Galaxy is made of the same dark matter that dominates the Galactic halo at large. In fact. 
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some years ago, we suggested that the supermassive compact dark object at the Galactic 
center may be a gravitationally stable ball of weakly interacting fermions in which the 
degeneracy pressure balances the gravitational attraction of the massive fermions (VioUier 
et al. 1992, 1993; VioUier 1994; Tsiklauri & VioUier 1996; BiUc, Munyaneza & VioUier 
1999). Such degenerate fermion balls (FBs) could have been formed in the early universe 
during a first-order gravitational phase transition (Bilic & VioUier 1997, 1998, 1999a,b). A 
further formation mechanism of FBs that is based on gravitational ejection of degenerate 
matter has recently been discussed in BiUc et al. 2000. 

There are three main reasons why it is worthwhile to study such degenerate FBs as an 
alternative to BHs at the center of the galaxies, in particular our own: 

(i) Introducing a weakly interacting fermion in the ~ 13 keV/c^ to ~ 17 keV/c^ mass 
range, one can explain the full range of the masses and radii of the supermassive 
compact dark objects, that have been observed so far at the galactic centers, in 
terms of degenerate FBs with masses ranging from 10^ to IO^-^Mq (Kormendy 
and Richstone 1995; Richstone et al. 1998). The maximal mass allowed for a FB 
composed of degenerate fermions of a given mass m/ and degeneracy factor gf is the 
Oppenheimer-Volkoff (OV) Umit Mqv = 0.54195 Mf>^ mj^ gj^ = 2.7821 x W^Mq 
(15 keV/mfC^)'^{2/gf)^ , where Mp£ = (hc/G)^ is the Planck mass (BiUc, Munyaneza 
& VioUier 1999). It is tempting to identify the mass of the most massive compact 
dark object ever observed at a center of a galaxy (Kormendy & Ho 2000), e.g. that of 
the center of M87, with the OV-Umit, i.e. Mqv = (3.2 ± 0.9) x lO^M© (Macchetto 
et al. 1997). This requires a fermion mass of 12.4 keV/c^ <^ rrif <^ 16.5 keV/c^ for gf 
= 2, or 10.4 keV/c2 ^ ruf ^ 13.9 keV/c^ for gf = 4. For Mqv = 3.2 x IO'^Mq such a 
relativistic FB would have a radius of Rqv — 4.45 Rqy — 1.36 mpc, where Rqy is the 
Schwarzschild radius of the mass Mqv- It would thus be virtually indistinguishable 
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from a BH, as the radius of the last stable orbit around a BH is 3 Rqy — 0.92 mpc 
anyway. The situation is quite different for a nonrelativistic FB of mass M = (2.6 ± 
0.2) X IO^Mq, which for the upper limit of the allowed fermion mass ranges, m/ = 
16.5 keV/c^ for gf — 2, or ruf — 13.9 keV/c^ for gf — 4, would have a radius bound 
by 16.7 mpc ^ R ^ 17.6 mpc, corresponding to ~ 7 x 10^ Schwarzschild radii, as the 
FB radius scales nonrelativistically like R oc ^f^^^ Sf"^^^ M~^/^. Such an object is far 
from being a black hole: its escape velocity from the center is ~ 1,700 km s~^. As 
the fermions interact only weakly with the baryons, baryonic stars could also move 
inside a FB without experiencing noticeable friction with the fermions (Tsiklauri and 
VioUier 1998a,b; Munyaneza, Tsiklauri and VioUier, 1998, 1999). Since the potential 
within ~ 10 mpc from the center is rather shallow, star formation in this region will 
be less inhibited by tidal forces than in the BH case. 

(ii) A FB with mass M = (2.6 ± 0.2) x lO^M© and radius R ^ 18.4 mpc is consistent 
with the current data on the proper motion of the stars in the central cluster around 
Sgr A*. This implies lower limits for the fermion masses of mj ^ 15.9 keV/c^ for g^ 

— 2 and m/ ^ 13.4 keV/c^ for g^ — 4, which partly overlap with the fermion mass 
ranges derived for M87. By increasing the fermion mass, one can interpolate between 
the FB and the BH scenarios. However, for fermion masses mj ^ 16.5 kcV/c^, for gf 

— 2 and m/ ^ 13.9 keV/c^ for gf — 4, the interpretation of some of the most massive 
compact dark objects in terms of degenerate FBs is no longer possible. It is quite 
remarkable that we can describe the two extreme cases, the supermassive compact 
dark object at the center of M87 and that of our Galaxy, in terms of self- gravitating 
degenerate FBs using a single fermion mass. This surprising fact is a consequence 
of the equation of state of degenerate fermionic matter; this would not be the case 
for degenerate bosonic matter. Indeed, for a supermassive object consisting of 
nonrelativistic self-gravitating degenerate bosons, mass and radius would scale, for a 
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constant boson mass, as R ex M~^, rather than R oc M~^/^, as for a supermassive 
object consisting of nonrelativistic self-gravitating degenerate fermions, for a constant 
fermion mass. The ratio of the radii of the supermassive objects with lO^'^M© and 
lO^-^M© would be 10^ in the boson case, instead of 10 as in the fermion case. Thus it 
would not be possible to fit mass and radius of both the supermassive compact dark 
object at the center of M87 and that of our Galaxy, in the boson case. We therefore 
conclude that, if we want to describe all the supermassive compact dark objects in 
terms of self-gravitating degenerate particles of the same kind and mass, these objects 
cannot be composed of bosons, they must consist of fermions. 

(iii) The FB scenario provides a natural cut-off of the emitted radiation at infrared 
frequencies 10^^ GHz, as is actually observed in the spectrum of the Galactic 
center (Bilic, Tsiklauri and Viollier 1998; Tsiklauri and Viollier 1999; Munyaneza and 
VioUier 1999). This is because matter, e.g. in the form of stars, gas, dust or dark 
matter, etc. falling from infinity at rest towards the FB, cannot acquire velocities 
larger than the escape velocity from the center of the FB, i.e. ~ 1,700 km s~^. 
Consequently, there is also a natural cut-off of the high-frequency tail of the radiation 
emitted by the accreted baryonic matter. This is quite a robust prediction of the FB 
scenario, because it is virtually independent of the details of the accretion model. In a 
thin disk accretion model, the radiation at the observed cut-off is emitted at distances 
~ 10 mpc from the center of the FB. This is also the region, where the gravitational 
potential becomes nearly harmonic due to the finite size of the FB. The fermion 
masses required for a cut-off at the observed frequency ~ 10^^ GHz depend somewhat 
on the accretion rate and the inclination angle of the disk assumed, but m/ ^ 20 
keV/c^ lor Qf — 2 ox ruf n keV/c^ for ^r^^ = 4 seem to be reasonable conservative 
upper limits (Tsiklauri and Viollier 1999, Munyaneza and Viollier 1999). 
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Summarizing the preceding arguments (i) to (iii), we can constrain the allowed fermion 
masses for the supermassive compact dark objects in our Galaxy to 15.9 keV/c^ ^ m/ ^ 
16.5 keV/c^ for gf = 2 or ISA keV/c^ ^ m/ ^ 13.9 keV/c^ for gf = 4, where the lower 
limits are determined from the proper motion of stars in the central cluster of our Galaxy, 
while the upper limits arise from the supermassive compact dark object at the center of 
M87. This fermion mass range is also consistent with the infrared cut-off of the radiation 
emitted by the accreted baryonic matter at the Galactic center. Of course, one of the major 
challenges will be to accommodate, within the FB scenario, the properties of Sgr A* which 
is perhaps peculiar to our galaxy. 

We now would like to identify a suitable candidate for the postulated weakly interacting 
fermion. This particle should have been either already observed, or its existence should 
have been at least predicted in recent elementary particle theories. The required fermion 
cannot be the gaugino-like neutralino, i.e. a linear combination of the bino, wino and the 
two higgsinos, as its mass is expected to be in the ~ 30 GeV/c^ to ~ 150 GeV/c^ range 
(Roszkowski 2001). It cannot be a standard neutrino either (however, see Giudice et al. 
2000), as this would violate the cosmological bound on neutrino mass and, more seriously, 
it would contradict the Superkamiokande data (Fukuda et al. 2000). However, the required 
fermion could be the sterile neutrino that has been recently suggested as a cold dark matter 
candidate in the mass range between ~ 1 keV/c^ to ~ 10 keV/c^ (Shi and Fuller 1999; 
Chun & Kim 1999; Tupper et al. 2000), although one would have to stretch the mass 
range a little bit and worry about the (possibly too rapid) radiative decay into a standard 
neutrino. This sterile neutrino is resonantly produced with a cold spectrum and near 
closure density, if the initial lepton asymmetry is ~ 10~^. Alternatively, it could be either 
the gravitino, postulated in supergravity theories with a mass in the ~ 1 keV/c^ to ~ 100 
GeV/c^ range (Lyth 1999), or the axino, with a mass in the range between ~ 10 keV/c^ and 
~ 100 keV/c^, as predicted by the supersymmetric extensions of the Peccei-Quinn solution 
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to the strong CP- problem (Goto & Yamaguchi 1992). In this scenario, the axino mass 
arises quite naturally as a radiative correction in a model with a no-scale supcrpotential. In 
summary, there are at least three promising candidates which have been recently predicted 
for completely different reasons in elementary particle theories. One of these particles could 
play the role of the weakly interacting fermion required for the supermassive compact dark 
objects at the centers of the galaxies and for cold or warm dark matter at large, if its mass 
is in the range between ~ 13 keV/c^ and ~ 17 keV/c^ and its contribution to the critical 
density is ~ 0.3. 



3. Outline of the paper 

The purpose of this paper is to compare the predictions of the BH and FB scenarios 
of the Galactic center, for the stars with the smallest projected distances to Sgr A*, based 
on the measurements of their positions during the last six years (Ghez et al. 2000). The 
projected orbits of three stars, SO-1 (SI), SO-2 (S2) and SO-4 (S4), show deviations from 
uniform motion on a straight line during the last six years, and they thus may contain 
nontrivial information about the potential. We do not rely on the accelerations determined 
directly from the data by Ghez et al. 2000, as this was done in the constant acceleration 
approximation which we think is not reliable. Indeed, the Newtonian predictions for the 
acceleration vary substantially, both in magnitude and direction, during the six years of 
observation. In view of this fact, we prefer to work with the raw data directly, trying to fit 
the projected positions in right ascension (RA) and dechnation of the stars in the BH and 
FB scenarios. For our analysis we have selected only two stars, SO-1 and SO-2, because their 
projected distances from SgrA* in 1995.53, 4.42 mpc and 5.83 mpc, respectively, make it 
most hkely that these could be orbiting within a FB of radius ~ 18 mpc. We thus may in 
principle distinguish between the BH and FB scenarios for these two stars. The third star. 
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SO-4, that had in 1995.53 a projected distance of 13.15 mpc from Sgr A*, and was moving 
away from Sgr A* at a projected velocity of ~ 990 km s~^, is now definitely outside a FB 
with a radius ~ 18 mpc. One would thus not be able to distinguish the two scenarios for a 
large part of the orbit of SO-4. 

In the following, we perform a detailed analysis of the orbits of the stars SO-1 and SO-2, 
based on the Ghez et al. 1998 and 2000 data, including the error bars of the measurements, 
and varying the unknown components of the position and velocity vectors of the stars in 
1995.4, z and v^. For simplicity, we assume throughout this paper that the supermassive 
compact dark object has a mass of 2.6 x IO^Mq, and is centered at the position of Sgr A* 
which is taken to be at a distance of 8 kpc from the sun. In fact, because of the small proper 
motion ^ 20 km s~^ of Sgr A*, there are strong dynamical reasons to assume in the BH 
scenario, that Sgr A* and the supermassive BH are at the same position, while in the FB 
scenario, Sgr A* and the FB could be off-center by a few mpc without affecting the results. 
We do not vary the mass of the supermassive compact dark object, as the calculations are 
not very sensitive to this parameter, as long as the mass is within the range of the error bar 
inferred from the statistical data on the proper motion of the stars in the central cluster 
(Ghez et al. 1998). 

This paper is organized as follows: In section 4, we present the main equations for the 
description of the supermassive compact dark object as a FB, as well as the formalism for 
the description of the dynamics of the stars in the gravitational field of a FB or a BH. We 
then investigate, in section 5, the dynamics of SO-1 and SO-2, based on the Ghez et al. 2000 
data, and conclude with a summary and outlook in section 6. 
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4. The dynamics of the steirs near the Galactic center 

As the stars near the Galactic center have projected velocities ^ 1,400 km s~^, one 
may very well describe their dynamics in terms of Newtonian mechanics for both the BH 
and the FB scenarios. Similarly, fermions of mass ruf ~ 13 keV/c^ to ~ 17 keV/c^, which 
are condensed in a degenerate FB of (2.6 ± 0.2) x lO^M©, are nonrelativistic, since their 
local Fermi velocity is certainly smaller than the escape velocity of ~ 1,700 km s~^ from the 
center of the FB. The fermions will, therefore, obey the equation of hydrostatic equilibrium, 
the Poisson equation and the nonrelativistic equation of state of degenerate fermionic 
matter 

Pf = Knf (1) 



with 



fc2 2\2/3 



X ^ __ _ 1 (2) 



Here, Pf and Uf, denote the local pressure and particle number density of the fermions, 
respectively. FBs have been discussed extensively in a number of papers (e.g. VioUier 1994; 
Bilic, Munyaneza and Viollier 1999; Tsiklauri and Viollier 1999). Here wc merely quote 
the equations that we need further below, in order to make this paper self-contained. The 
gravitational potential of a degenerate FB is given by 



$(r) = < 



GM 

, X > Xq , 



ax 

where a is an appropriate unit of length 

Snh^ /l5keV\^/^ 

mj gf Mq' 

r = ax is the distance from the center of the FB and R — axo the radius of the FB. The 
dimensionless quantity v{x), that is related to the gravitational potential $(r) through 
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eq.(3), obeys the Lane-Emden differential equation 



(5) 



dx^ 2.1/2 ' 

with poly tropic index n — 3/2. For a pure FB without a gravitational point source at the 
center, the boundary conditions at the center and the surface of the FB are v{0) — v{xo) — 0. 
All the relevant quantities of the FB can be expressed in terms of v and x, e.g. the matter 
density as 

where m/ and gf are the mass and the spin degeneracy factor of the fermions and 
antifermions, respectively, i.e. gf — 2 for Majorana and gf — 4 for Dirac fermions and 
antifermions. Based on eqs.(5) and (6), the mass enclosed within a radius r in a FB is given 

by 

M(r) = / Aw p r'^ dr = - Mq (v'(x)x - v(x)) , (7) 
and the total mass of the FB by 

M = M(R) = - M© v'{xo) xo . (8) 

From eq.(5), one can derive a scaling relation for the mass and radius of a nonrelativistic 
FB, i.e. 



MR^ = xo\v'{xo)\xla^ Mq 



91.869 /1^ f2V 



G^rrA \gf, 



Here v{x) is the solution of eq.(5) with v(0) = and v'(0) = 1, yielding v{xq) — again at 
xq — 3.65375, and v'{xo) — - 0.742813. The precise index of the power law of the scaling 
relationship (9) depends on the polytropic index of the equation of state (1). As the mass of 
the FB approaches the OV limit, this scaling law is no longer valid, because the degenerate 
fermion gas has to be described by the correct relativistic equation of state and Einstein's 
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equations for the gravitational field and hydrostatic equilibrium (Bilic, Munyaneza & 
VioUier 1999). 

We now turn to the description of the dynamics of the stars near the Galactic center. 
The mass of the BH and FB is taken to be M = 2.6 x lO^M©. In order to emphasize 
the differences between the FB and the BH scenarios, we choose the fermion masses ruf — 
15.92 keV/c^ ior gf — 2 or m,f — 13.39 keV/c^ for gf — 4. These are the minimal fermion 
masses consistent with the mass distribution inferred from the statistics of proper motions 
of the stars in the central cluster (Munyaneza, Tsiklauri and VioUier, 1999; Ghez et al. 
1998). The dynamics of the stars in the gravitational field of the supermassive compact 
dark object can be calculated solving Newton's equations of motion 

r = ^ r , (10) 

taking into account the position and velocity vectors at e.g. to — 1995.4 yr, i.e. 
f{to) = {x,y,z) and f(to) = {vx-,Vy,Vz). For the FB scenario, M{r) is given by eq.(7), while 
in the BH case it is replaced by M of eq.(8). The ,x-axis is chosen in the direction opposite 
to the right ascension (RA), the y-axis in the direction of the declination, and the 2;-axis 
points towards the sun. The BH and the center of the FB are assumed to be at the position 
of Sgr A* which is also the origin of the coordinate system at an assumed distance of 8 kpc 
from the sun. 



5. Analysis of the orbits of SO-1 and SO-2 

In 1995.4, the projected positions and velocities of SO-1 reported by Ghez et al. 1998, 
were x = -0.107", y = 0.039", = (470 ± 130) km s'^ and Vy = (-1330 ± 140) km s-^ 
We now investigate how the projected orbits, calculated using eq.(lO), are affected by (i) 
the error bars of Vx and Vy of SO-1 measured in 1995.4, (ii) the lack of knowledge of z of 
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SO-1 in 1995.4, (iii) the lack of information on of SO-1 in 1995.4. We then compare the 
results with the SO-1 data recently reported by Ghez et al. 2000. 

Fig.l shows the RA of SO-1 as a function of time, taking into account the error bars 
of Vx and Vy and choosing z = = in 1995.4. The top panel represents the RA of SO-1 
in the BH scenario, while the bottom panel illustrates the same quantities in the FB case. 
Prom Fig.l we conclude that, ior z — — in 1995.4, the error bars of Vx and Vy of 1995.4 
do not allow for a fit of the new Ghez et al. 2000 data of SO-1 in the BH scenario, whereas 
the data are described quite easily within the error bars in the FB case. In Fig. 2, the 
declination is plotted as a function of time for the same values of Vx, Vy, z and Vz in 1995.4. 
We arrive at the same conclusion as in Fig.l: For z — Vz — the error bars of Vx and Vy 
allow for a fit of the data in the FB scenario only. 

As a next step, the sensitivity of the orbits to the z- coordinate of SO-1 in 1995.4 is 
investigated. To this end, we restrict ourselves to bound orbits of SO-1 only. The conserved 
total energy of the star SO-1 is given by 

1 l2 

E ^ -mr + m$(r) , (11) 

where the unknown star mass can be chosen as m = 1 without loss of generality. SO-1 is 
unlikely to have a total energy £^ > 0, because, in the absence of swing-by events caused 
by stars of the central cluster, SO-1 will have to fall in with a velocity that is inconsistent 
with the velocity dispersion of the stars at infinity. The condition < thus yields upper 
limits, li'^l < \^^\ and j^l < l^ool, which depend on Vx and Vy as can be seen from Table 1. 
In this context, it is worthwhile to note that, at a radius equal to the projected distance of 
SO-1 to Sgr A* in 1995.4, the escape velocity from a BH is 2, 250 km s~^, while that from a 
FB is 1, 613 km s~^. The escape velocity from the center of the FB is 1,672 km s~^. 

Fig. 3 presents the sensitivity of the RA of SO-1 with respect to the z-coordinate in 
both the BH and FB scenarios. In the case of a BH, the RA depends strongly on the value 
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of z in 1995.4, while the 2;-dependence in the FB scenario is rather weak. In both the top 
and bottom panels, = 340 km s'^, Vy = —1190 kms~^ and = has been assumed 
while z is varied, all in 1995.4 values. In the BH scenario, none of the other input values 
for Vx and Vy would fit the new Ghez et al. 2000 data if we restrict ourselves to bound 
orbits. In the FB case, the input values — 340 km and Vy — —1470 km s~^ describe 
the Ghez et al. 2000 data as well as = 340 km s~^ and Vy = —1190 km s^^. Looking 
at Fig. 3, we thus conclude that the RA data of SO-1 are well fitted with \z\ ~ 0.25" in 
the BH scenario, and with \z\ ^ 0.1" in the FB case. The dependence of the declination on 
z is shown in Fig.4. In order to describe the Ghez et al. 2000 data for the declination of 
SO-1, we require 0.25" ^ \z\ ^ \zoo\ = 0.359" for the BH scenario, while in the FB case, the 
declination can be fitted with \z\ ^ 0.359". 

Fig.5 represents an investigation of how the RA of SO-1 is affected by of 1995.4. 
In this graph we have chosen Vx — 340 km s~^, Vy — - 1190 km and z — in 1995.4. 
Increasing \vz\ up to its maximal value = 1879 km does not help fitting the RA data 
of SO-1 in the BH scenario. In the FB case, the orbits are rather insensitive to If^l- Thus 
all \vz\ ^ \v'^\ — 1036 km s~^ fit the RA data of SO-1 quite well. The weak dependence 
of the RA on l^;^;! in the FB scenario is due to the harmonic oscillator hke shape of the 
FB potential at small distances, where the Newtonian equations of motion nearly decouple 
in Cartesian coordinates. Fig. 6 exhibits the declination as a function of time for various 
values of \vz\, keeping Vx, Vy and z in 1995.4 as in Fig.5. The top panel of Fig. 6 shows that 
increasing l^^l from zero up to its maximal value = 1879 km s~^ barely helps fitting 
the data in the BH scenario. In the FB scenario, the declination may be described by \vz\ ^ 
900 km s-^ 

Summarizing the results of Figs. 3-6, we can plot in Fig.7 the z — phase-space of 
1995.4 that fits the data. The small range of acceptable \z\ and \vz\ values in the BH 
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scenario (solid vertical line) reflects the fact that the orbits of SO-1 depend strongly on 
these two parameters. Conversely, the weak dependence of the orbits on and \vz\ in the 
FB case is the reason for the much larger z — phase-space that is fitting the Ghez et al. 
2000 data of SO-1, as shown by the dashed box. The dashed and sohd curves describe the 
E — (just bound) orbit in the FB and BH scenarios, respectively. 

Fig.8 shows some typical projected orbits of SO-1 in the BH and FB scenarios. The 
Ghez et al. 2000 data of SO-1 may be fitted in both scenarios with appropriate choices of v^, 
Vy, z and in 1995.4. The inclination angles of the orbit's plane 6 = arccos (^Lz/\L\^, with 

— * 

L — mf X r, are shown next to the orbits. The minimal inclination angle that describes 
the data in the BH case is ^ = 70°, while in the FB scenario it is ^ = 0°. In the BH 
case, the minimal and maximal distances from Sgr A* are = 0.25" and r^ax = 0.77", 
respectively, for the orbit with z = 0.25" and = which has a period oiT ^ 161 yr. The 
orbits with z — 0.25" and Vz — 400 km s~^ or z — 0.25" and Vz — 700 km s""^ have periods 
of T fa 268 yr or T 3291 yr, respectively. In the FB scenario, the open orbit with z — 0.1" 
and f z = has a "period" of T ^ 77 yr with rmin = 0.13" and Vmax = 0-56". The open 
orbits with z — 0.1" and — 400 km s~^ or z — 0.1" and — 900 km s~^ have "periods" 
oiT Ki 100 yr or T fs 1436 yr, respectively. 

Fig.9 is a prediction of \vz\ as a function of time for both the BH and FB scenarios and 
various acceptable input parameters (see Fig. 7). This shows that the radial velocity \vz\ of 
SO-1, if measured in a few years time, could serve to distinguish between the BH and the FB 
scenarios. For the BH case, we predict by the year 2005 a radial velocity l^^l ^ 900 km s~^, 
while in the FB scenario the radial velocity will be l^^l ^ 500 km s~^. Radial velocities 
\vz\ ^ 1000 km s~^ before 2010 would be excluded in both the BH and FB scenarios. 

We now repeat the analysis in the case of the star SO-2. The x- and y-components 
of the position and velocity vectors of SO-2 at to = 1995.4 yr are f{to) = (0.007", 0.151") 
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and f{to) = (-290 ± 110, -500 ± 50) km s-^ respectively (Ghez et al. 1998). Restricting 
ourselves to bound orbits, the Ghez et al. 2000 data of SO-2 can be fitted in the BH 
scenario, with Vx = -290 km s~^, Vy = -500 km s~^, \z\ ~ 0.25" and \vz\ ^ 1280 km s~^ in 
1995.4. In the FB scenario, the allowed z — phase-space is ^0.1" and \vz\ ^ 1000 
km s~^ for the same values of and Vy in 1995.4. The range of acceptable values oi \z\ and 
1^2! is shown in Fig. 10, where the solid and dashed curves denote the limits on \z\ and l^^l 
for £^ = (just bound) orbits in the BH and FB cases, respectively. Here again the z — Vz 
phase-space turns out to be much larger in the FB (dashed box) than in the BH scenario 
(vertical sohd line). 

Fig. 11 exhibits some typical projected orbits of SO-2 corresponding to the acceptable 
z — Vz phase-space of Fig. 10 in both the BH and FB scenarios. Next to the orbits the 
inclination angles of the orbit's plane are shown. In the BH scenario, the lower limit of the 
inclination angle is ^ = 74°, for the orbit (1) with z — 0.25" and Vz — 0. The minimal and 
maximal distances from the center are rmin — 0.045" and Vmax — 0.307", respectively. The 
period of this closed orbit is minimal with T ^ 34 yr. However, this value disagrees with 
the minimal period T !^ 16 yr quoted by Ghez et al. 2000. The orbit (2) with z — 0.25" and 
Vz — 500 km s~^ or 2; = 0.25" and Vz =1000 km s~^ have periods T 42 yr or T 135 yr, 
respectively. In the FB scenario the minimal inclination angle of the orbit's plane is ^ = 0°. 
The open orbit (1) with z = 0.25" and Vz = has minimal and maximal distances from the 
center of rmin — 0.15" and Vmax — 0.31", respectively, with a "period" of T 46 yr. The 
open orbits with (1) z = 0.25" and Vz = 500 km s"^ or (2) z = 0.25" and Vz = 1000 km s"^ 
have "periods" T pa 51 yr or T pa 174 yr, respectively. The minimal "period" T pa 37 yr is 
obtained for an open orbit with z — Vz — and inclination angle ^ = 0°. 

In Fig. 12 the predictions for the radial velocity \vz\ are plotted as a function of time 
for various acceptable input values of z and Vz in 1995.4. It is seen that a radial velocity of 
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\vz\ ^ 1000 km s ^ is excluded in the FB scenario. 

6. Summciry and outlook 

In this paper, we have shown that the orbits of SO-1 and SO-2, observed by Ghez et al. 
2000 during the last six years, arc consistent with cither a BH or a FB of 2.6 x lO^M© at 
the center of the Galaxy. In order to fit the data in the BH scenario, SO-1 and SO-2 must 
have had in 1995.4 a z-coordinate of \z\ 0.25" and radial velocities of \vz\ ^ 750 km s~^ 
oi' 1^2 1 ^ 1280 km s~^ for SO-1 or SO-2, respectively In the BH scenario, the orbits of SO-1 
and SO-2 strongly depend on z in 1995.4. The new data of SO-1 and SO-2 can be fitted in the 
FB scenario with z — — in 1995.4. Due to the weak dependence of the orbits on and 
\vz\ in the FB case, \z\ ^ 0.1" and {v^l ^ 900 km s"^ for SO-1 or \z\ ^ 0.25" and \vz\ ^ 1000 
km s~^ for SO-2 are also consistent with the Ghez et al. 2000 data. As new measurements 
become available, the acceptable z — Vz phase-space of 1995.4 could be further constrained. 

We have plotted some typical orbits of SO-1 and SO-2 in both the BH and FB scenarios 
and have shown that the minimal inclination angle of the orbit's plane is ^ = 70° as in 
the BH case and ^ = 0° in the FB scenario. We have established that by the year 2005, 
the measurement of the radial velocities \vz\ of both SO-1 and SO-2 could discriminate 
between the two scenarios of the supermassive compact dark object at the Galactic center. 
In concluding, it is important to note again that, based on the Ghez et al. 2000 data of 
the stars SO-1 and SO-2, the FB scenario cannot be ruled out. On the contrary, in view of 
the z — Vz phase-space, that is much larger in FB scenario than in the BH case, there is 
reason to treat the FB scenario of the supermassive compact dark object at the center of 
our Galaxy with the respect it deserves. 

We now turn to the discussion of promising techniques for proving or disproving the 
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FB or BH scenarios of the Galactic center. These can be basically divided into two classes: 

(i) probing the gravitational potential at distances ^ 0.5" from the Galactic center, and 

(ii) observing the decay or the annihilation of the fermions of the FB into visible particles, 
at distances ^ 0.5" from the Galactic center. 

In the first category, the most promising method is still monitoring the proper motions 
and radial velocities of stars that are located at projected distances ^ 0.5" from Sgr A*, 
as well as interpreting these observations in terms of the FB or BH scenarios. However, a 
further interesting possibility is the observation of gravitationally lensed stars in the line 
of sight behind the FB (Bihc, Nikolic and VioUier 2000), as a FB of ~ 18 mpc radius is a 
much more efficient gravitational lens than a BH of the same mass. In fact, a star crossing 
the line of sight with a minimal distance of ~ 0.2 mpc, ~ 200 pc behind the center of the 
FB, will produce for a few years up to three distinct moving images within or just outside 
the Einstein ring radius of ~ 0.13" or ~ 5 mpc. Two of these star images will pop out of 
nothing at some point inside the Einstein ring, and they will first separate, then approach 
each other again and finally annihilate each other at a difi'crcnt point within the Einstein 
ring. The third image moves around the Einstein ring while the projection of the source 
crosses the ring area. Unfortunately, the rate for such a remarkable event, which would be 
the smoking gun for an extended supermassive object, is estimated to be only ~ 10~^/yr. 
In the BH case, the manifestations of gravitational lensing are less spectacular with only 
two observable lensed moving images (Wardle and Yusef-Zadeh 1992). Another possible 
test of the gravitational potential could be the spectrum of the radiation emitted by the 
accreted baryonic matter, once the model dependence of the calculations can be controlled. 

In the second category, the particle content of FBs could be proven e.g. through the 
radiative decay of the fermion (assumed here to be a sterile neutrino) into a standard 
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neutrino, i.e. / — > v^. If the lifetime for this decay is 1.4 x 10^^ yr, the luminosity of the 
FB would be 4 X 10^^ erg/sec. The signal would be a sharp X-ray hne of ~ 8 keV for gf 
= 2 or ~ 7 keV for gf = 4. The X-ray luminosity would be tracing the fermion matter 
density. Of course the spatial resolution of the X- ray telescope would have to be ^ 0.5". 
A further possible test could be the annihilation of two fermions into two or three photons, 
i.e. // — >■ 77 or // — > 777. However, the branching ratios of these two annihilation 
processes with respect to a presumably dominant but unobservable // v'v' channel are 
most probably too small to be observable. Nevertheless, the signal would be a sharp line of 
~ 16 keV ((7/ = 2) of ~ 14 keV {gf = 4) for // ^pj, and a continuous spectrum with a 
maximum at ~ 10 keV [gf = 2) or ~ 9 keV [gf = 4) for // 777- For s-wave annihilation 
the X-ray luminosity would trace the square of the fermion matter density, while for p-wave 
annihilation the concentration of the X-ray emission towards the center of the FB would be 
even more pronounced. 
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Figure captions: 

Fig.l : The RA of SO-1 as a function of time, taking into account the error bars of the 
projected velocity in 1995.4 (Ghez et al. 98). The top panel is drawn for the case of a black 
hole (BH), while the bottom panel corresponds to a fermion ball (FB), both centered at 
Sgr A*. The labels of the graphs stand for the velocity components: (1) — 470 km s~^, 
Vy — -1330 km s"-*^ (median values); (2) Vx — 340 km s"-*^, Vy — -1190 km s"-*^; (3) Vx — 340 
km s-^ Vy = -1470 km s^^ (4) = 600 km s-^ Vy = -1190 km s^^; (5) Vx = 600 km s-\ 
Vy — -1470 km s~^, all for 1995.4. In this plot, z and v^ are taken to be zero in 1995.4. One 
thus concludes that, within the error bars of the projected velocities of 1995.4, the data 
points recently reported by Ghez et al. 2000 can be fitted in the FB scenario with v^ = z = 
(lines 2 and 3), while this is not the case in the BH scenario. In this figure, the masses of 
both the FB and the BH, are M = 2.6 x lO^M©. The fermion mass is taken to be m/c^ 
= 15.92 keV with a degeneracy factor of gf — 2, or m/c^ = 13.39 keV with gf — 4. Fig. 2: 
This is an analysis similar to that of Fig.l but for the declination of SO-1 as a function 
of time. Assuming v^ — z — in 1995.4, the data points can be represented within the 
error bars of Vx and Vy in 1995.4 in the FB scenario only. Fig. 3: Here we investigate the 
dependence of the RA of SO-1 on the 2;-coordinate. The data can be described with Vx = 
340 km s~^, Vy = -1190 km s~^, v^ = and \z\ ~ 0.25" in the BH case. In the FB scenario, 
the data can be fitted for the same values of Vx, Vy, Vz and \z\ ^ 0.1". We note that the 
RA strongly depends on z in the BH case, while in the FB scenario the dependence on z 
is rather weak, as the equations of motion nearly decouple for small x, y, and z. The limit 
value l^ool = 0.359" corresponds to a just bound orbit with total energy E — and Vz — 
(see Table 1). Fig.4 : The dependence of the declination of SO-1 on \z\ for Vz — 0, Vx — 
340 km s~^ and Vy = -1190 km s~^, all in 1995.4 values. The declination can be fitted from 
\z\ a; 0.25" to \zoo\ = 0.359" in the BH case. In the FB scenario all values \z\ ^ 0.359" fit the 
data. Fig. 5: Here we explore how the RA of SO-1 depends on l^^l while keeping z — 0, Vx — 
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340 km and Vy — -1190 km s~^, taken in 1995.4. Increasing \vz\ up to its maximal value 
= 1879 km s"-*^, which corresponds to a £^ = (just bound) orbit with z = 0, does not 
help fitting the data in the BH scenario. In the FB case, there is a weak dependence on v^, 
and all the values \vz\ ^ 1879 km describe the data. Fig. 6: The declination of SO-1 as 
a function of \vz\ for z — — 340 km s~^ and Vy — -1190 km s~^, all in 1995.4 values. 
In the BH case, the data can only be described for very large l^^l, while in the FB scenario 
the data are fitted for l^^l ^ 900 km s~^. Fig. 7: The z — phase-space of 1995.4 that fits 
the SO-1 data with — 340 km s~^ and Vy — - 1190 km s~^ in 1995.4. In the BH case,the 
data require \z\ ~ 0.25" (vertical solid line), as the orbits strongly depend on z, and \vz\ ^ 
750 km s~^. The large range of acceptable \z\ and \vz\ values in the FB case, denoted by 
the dashed rectangular box, is due to the fact that the orbits depend weakly on these two 
parameters. Values of \z\ ^ 0.1" and l^^j ^ 900 km s~^ are consistent with the measured 
positions of SO-1 in the FB scenario. The solid and dashed curves correspond to the limits 
of possible values of and for £^ = (just bound) orbits of SO-1 in the BH and FB 
cases, respectively. Fig. 8: Examples of typical orbits of SO-1 for — 340 km s"^ and Vy — 
- 1190 km s~^ in 1995.4. The inclination angles of the plane of the orbit are indicated next 
to the curves. In the BH case, the orbit with z = 0.25" and = in 1995.4 has minimal 
and maximal distances from Sgr A* of rmin — 0.25" and Vmax — 0.77" , respectively, with a 
period of T 161 yr. The orbits with z — 0.25" and — 400 km s~^ ot z — 0.25" and 
= 700 km s~^ in 1995.4 have periods of T ^ 268 yr or T ^ 3291 yr, respectively. In the 
FB scenario, the open orbit with z — 0.1" and Vz — has a "period" of T w 77 yr with 
r^.„ = 0.13" and Vmax = 0.56". The open orbit with z = 0.1" and Vz = 400 km s"^ in 1995.4 
has a "period" T fa 100 yr, while that with z = 0.1" and Vz = 900 km s'^ has a "period" 
of T fti 1436 yr. We note that the minimal inclination angle of the orbit's plane is ^ = 70° 
in the BH case and ^ = 0° in the FB scenario. Fig. 9: The prediction for \vz\ as a function 
of time for SO-1. In the BH scenario with l^j = 0.25" and i^;;; = in 1995.4, \vz\ should be 
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^ 900 km by the year 2005. The FB scenario predicts a \vz\ ^ 500 km by the year 
2005, for Vz = and \z\ = 0.1" in 1995.4. Thus, once the radial velocities are measured, 
this figure could serve to distinguish between the BH and FB scenarios of Sgr A*. Fig. 10: 
The z — Vz phase-space of 1995.4 that fits the SO-2 data with Vx — - 290 km s~^ and Vy — 
- 500 km s~^ in 1995.4. The sohd and dashed curves denote the hmits on \z\ and \vz\ for 
E = (just bound) orbits in the BH and FB scenarios, respectively. In the BH case, the 
vertical solid line at \z\ ~ 0.25" up the BH curve stands for the region of the allowed values 
of 1 2; I and \vz\ for SO-2. In the FB scenario, the large allowed z — Vz phase-space denoted by 
the dashed rectangular box, is due to the fact that the open orbits depend rather weakly 
on these two parameters. Fig.ll: Examples of typical orbits of SO-2 for = - 290 km s~^ 
and Vy — - 500 km s~^ in 1995.4. The labels of the orbits in both the upper (BH) and lower 
(FB) panels stand for: (1) z = 0.25", Vz = 0; (2) z = 0.25", Vz = 500 km s"^; (3) z = 0.25", 
Vz = 1000 km s~^. The parameters of the label (0) z = 0, Vz = 150 km s^^ fit the SO-2 data 
in the FB scenario only. The inclination angles of the plane of the orbit are indicated next 
to the curves. In the FB case, these values correspond to the maximal possible \z\ and \vz\ 
in 1995.4 for bound orbits. Bearing in mind that \z\ ^ 0.25" in the FB case, the inclination 
angle can take any value from 6* = 0° to 6' = 82°. As an example, an orbit with 9 = 29° 
with the label (0) is shown in the lower panel. However, in the BH scenario, the lower limit 
for the inchnation angle is ^ = 74°. The closed orbit with label (1) has r^m — 0.045", 
Tjnax = 0.307" and a period of T ^ 34 yr. The closed orbits with labels (2) and (3) have 
periods of T 42 yr and T w 135 yr, respectively. The open orbit with label (1) has a 
"period" of T 46 yr with rmin — 0.15" and Vmax — 0.31". The open orbits with labels (2) 
and (3) have "periods" of T fa 51 yr and T fa 174 yr, respectively. Fig. 12: The prediction 
for 1^2! as a function of time for SO-2. This could serve to distinguish between the FB and 
the BH scenarios of the supermassive dark object at the Galactic center. A radial velocity 
Vz ^ 1000 km would rule out the FB scenario. 
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